Non-standard interaction in neutrino oscillations and recent Daya Bay, T2K 

experiments 



Rathin Adhikari,^^0 Sabyasachi Chakraborty,'^'[l| Arnab Dasgupta,^'[l and Sourov Roy^'[§ 

Centre for Theoretical Physics, Jamia Millia Islamta (Central University), Jamia Nagar, New Delhi-110025, INDIA 
'^Department of Theoretical Physics, Indian Association for the Cultivation of Science, 
2 A & 2B Raja S.C.Mulhck Road, Jadavpur, Kolkata 700 032, INDIA 

We study the possible constraints on non-standard interaction(NSIs) in a model independent way 
by considering the recent results from T2K and Daya Bay neutrino oscillations experiments. Using 
perturbation method we present generic formulas (suitable for T2K baseline and for large 6^13 as 
evident from Daya Bay) for the probability of oscillation for i/^ Ve, taking into account NSIs at 
source (e"), detector (e"*) and during propagation (e"*) of neutrinos through matter. Two separate 
cases of perturbation with small (slightly large) NSI (e™^ ~ 0.03(0.18)) are discussed in detail. Using 
various possible presently allowed NSI values we reanalyze numerically the S13 — 5 allowed region 
given by recent T2K experimental data. We obtain model independent constraints on NSIs in the 
5 — e™^ plane using the ^13 value as measured by Daya Bay, where 5 is the CP violating phase. 
Depending on 5 values significant constraints on ter and Ett, in particular, are possible for both 
hierarchies of neutrino masses. Corresponding to T2K's 66% confidence level result, the constraints 
on Err is shown to be independent of any 5 value. 

PACS numbers: 13.15.-|-g, 14.60.Pq, 14.60.St 



I. INTRODUCTION 

Neutrino oscillations successfully describe neutrino 
flavor transitions. The recent superbeam and reactor 
neutrino experiments have provided enormous insights 
to unravel the exact value of the vacuum mixing angle 
6'i3. To emphasize this point, the T2K[i| experiment 
observed indications of i'^ — )■ Ve appearance, by pro- 
ducing a conventional neutrino beam at J-PARC and 
directed 2.5° off- axis to a detector situating at 295 Km 
away. The bounds on ^13 which, T2K came up with was 
0.03(0.04) < sin^ 26I13 < 0.28(0.34) for (5 = and normal 
(inverted) hierarchy. The reactor neutrino experiments 
like Daya BayQ and RenoQ provided compelling 
evidences for a relatively large angle ^13, with 5.2(t and 
4.9tT results respectively. These recent reactor neutrino 
results indicate Oiz very close to 8.8°. 

In this work we considered non-standard interactions 
(NSIs), occurring from four-fermion operators. In addi- 
tion to the standard model Lagrangian density, we con- 
sider the following non-standard interactions in the low 
energy effective theory during the propagation of neutri- 
nos through matter 

C'j^Si - -2V2Gj.e^^ [h'-Pf] [i^cI^Pl^p] , (1) 

where f = e,u,d and P — Pl,Pr where Pz, = (1 — 75)/2 
and Pr, = (1 -I- 75)72. In our subsequent sections these 
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non-standard interactions (NSIs) are related to e™^ as 
mentioned later. In the neutrino oscillation experiments 
for the NSIs at the source and detector the following 
lagrangian densities as low energy effective theory corre- 
sponding to charged current interactions due to leptons 
and quarks may be considered: 

Cf^^j - -2V2GFe^?^ [hrPla] [>>^7^.Pl'^s] (2) 

Cfj§j = -2V2GFe"rKd [u-ff^Pd] [l^j^^PLiys] 

+h.c (3) 

In our subsequent sections both these NSIs - e^^''^ 
and e"^^ contribute to e^^ and e[^^ corresponding to 
the appropriate interactions at source and detector for 
the neutrinos respectively. From the Lagrangian, we 
observe that the NSI parameters do not possess any 
mass dimension. However, if NSIs are related to the 
underlying new physics, then they should be considered 
as a first order term in the perturbation series and not in 
the zeroth order [4|. To reiterate this, we know the NSI 
parameters are related to the new physics scale in the 
form e ^ [Mw /M^qsi)'^ , where M^si signifies the new 
physics scale. So if we consider the new physics scale to 
be around a few TeV, then the NSI parameters should 
not be greater than a few percent. In general the NSI 
parameters can be categorized into two different parts. 
One is the NSI during propagation, and the other being 
NSI at the source and at the detector. It is worthwhile 
to note that the present bounds on the NSI parameters 
during propagation are not very stringent 

Non Standard interaction and its implications in a 
model independent way, as well as in different models 
have already been studied very extensively in the 
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literature [6l-[l8|. Many authors have studied their 
impact on solar neutrinos [l9l - l22j . atmospheric neu- 
trinos f20 42q | . conventional and upgraded neutrino 
beamsMna^, [H, H, HO], neutrino factories 0, i, [13- 
[l^. [l^lsiMSSj l ■ beta beams [341] . supernova neutrinos 
|35l - l37| . cosmological relic neutrinos [Sl] and neutrino- 
instability problem [39}. Also in some cases similarity 
of such effective interactions with CPT violations [40] 
and probing such interactions at LHC [4l| have been 
discussed. In the context of solar neutrinos possible 
confusion of nonzero mixing angle ^13 in the presence of 



NSI and hint of NSI have been mentioned [42|. One of 
the most striking features of NSI parameters, is to cloud 
the sensitivity of ^13 by orders of magnitude, which was 
shown very explicitly in [isl l43j . for neutrino factories 
and reactor neutrino experiments respectively. 

To elaborate the plan of our paper, at first in section 
II, we present a generalized prescription (suitable for rel- 
atively short baseline of T2K) by following the works of 
pil - liq . which in the literature is also known as the '-y/e 

method of perturbation theory', where e = ^ 0.03 

and represent a mathematical formulation by consider- 
ing a relatively large sin^is ~ ~ 0.18. We divide the 
Hamiltonian, consisting of the standard matter interac- 
tion and NSI during propagation, into zeroth order part 
and a perturbative part, where \fi is the perturbation 
parameter. Our next task is to compute the S-niatrix 
elements from these Hamiltonians. In section III, we in- 
voke the idea of NSI parameters at the source and the 
detector. Previous bounds on these parameters were con- 
strained by lepton and pion decays [13, [13], which were 
of the order of 0(0.1). However, the present bounds on 
NSI parameters at the source and at the detector are very 
strong[3]. Due to this reason we assume the NSI param- 
eters, present at the source and the detector are of the 
order of e. In section IV and V we have considered two dif- 
ferent cases, one with the consideration of e 



Q/3 



\fk and 

the other case with e™^ ~ e, where e™^ is the NSI param- 
eter during propagation. We have presented the expres- 
sion of the probability up to second order in e, by taking 
into consideration all these effects, such as the standard 
matter interaction, NSI during propagation and NSI at 
the source and at the detector. We were able to match 
the results obtained from the analytical expressions with 
that of the full numerical study. This also shows the re- 
markable power of this perturbation method. We show 
that because of the presence of NSIs at the source and at 
the detector, one can have a non-zero oscillation proba- 
bility at the source itself without the neutrino traversing 



any length. This is coined as the zero distance effect [13- 
[SOf or the near detector effect [Slj in the literature. This 
effect is a manifestation of the non-unitarity of the mix- 
ing matrix by considering NSIs at the source and at the 
detector. 

It is also important to note that, in principle one can 
also follow the method of matrix perturbation, to obtain 
the expression for the probability. In that case one has 
to compute the modification of the PMNS matrix[53l.[53|. 
due to the inclusion of the standard matter interaction 
and NSI during propagation. The modified PMNS ma- 
trix has to be diagonalised, the eigenvectors and the 
eigenvalues are to be extracted from the modified PMNS 
matrix. After that NSI at the source and at the detec- 
tor are to be included, to compute the overall expression 
of the oscillation probability. Similar approaches were 
followed by the authors of [s^, H^l 

In section VI, using Daya Bay and T2K experimental 
results we have discussed numerical analysis in obtain- 
ing the constraints in the 5-NSI plane. Here, the larger 
model-independent allowed values of NSI (not considered 
in our perturbative approach) have been considered for 
the analysis. 



II. MATHEMATICAL FORMULATION FOR 
LARGE 6I13 PERTURBATION THEORY 

The recent reactor based neutrino experiments have 
provided substantial proof for a relatively larger ^13. 
Based on the works of [43 - |46| . we describe a mathemat- 
ical prescription to show the effects of the non-standard 
interactions during propagation in neutrino oscillations. 
We consider the channel — > z^e, as followed by the 
recently concluded T2K experiment. Using the present 
experimental values of 6113 and the mass squared differ- 
ences, we formulate 



sm6'i3 = si3 



'21 



0.03. 



(4) 



This section elaborates the basic principles of our pertur- 
bative approach. In the Shrodinger picture, a neutrino 
with flavor a obeys the evolution equation [57| 



Km 



(5) 



where the Hamiltonian (after extracting constant diago- 
nal matrix irrelevant for flavor transition as it generates 
a phase common to all flavors) is given as 
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The inclusion of the standard matter effect to the Hamil- 
tonian is commonly known as the MSW effect (ssl [59| . 
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is the PMNS[53l |5J] matrix in vacuum. A = 2EVcc 
represents the interaction of the neutrino with matter, 
more precisely with electrons. E is the energy of the 
neutrino, Vcc represents the charge current interaction 
and given by Vcc = V^GpNc, where Gp is the Fermi 
coupling constant and N,, is the electron number density. 
By taking Am^i outside the square brackets, from Eq. 



(j6|), we redefine the matter interaction as A 
and define a 
A = 0.06 ~ e. 
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A/Aml^ 

e. For the T2K experiment. 



becomes. 
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where we have defined J723 in Eq.Q- This Hamiltonian 
in the tilde basis can now be written as a sum of the 
Hamiltonians of different orders {T-Lm = Hq + Hi), where 
the ordering is done with respect to ^/e. For example 
the zeroth order Hamiltonian, as a function of standard 
matter interaction looks like. 



From Eq.®, let us first consider the case where NSI 
is absent. As a method to simplify calculations, it is 
convenient to work in the tilde basis, which we define 

as i^a = (C^23)a/3^/3- ^^^^ basis the Hamiltonian, con- 
sisting only the standard matter interaction part, or Hm 



Ho 



and similarly. 
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Here Hi is the perturbed part of the Hamiltonian in 
standard matter. The different matrices in the perturbed 
Hamiltonian in the tilde basis comprises of four differ- 
ent orders in y^, which are -y/e, e, , and respectively. 

Now we include the NSI matrix during propagation. 
The Hamiltonian consisting of these NSI parameters has 
the formH, 
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where, 



E 



jpni 

-a/3 



(12) 



where nj is the number density of the fermion / Q . Here 
e™^ , {a,l3 — e, /i, r) are non-standard interaction (NSI) 
parameters of neutrinos, propagating through matter, de- 



fined as e 



al3 



To include the Ti.NSi matrix 
in the perturbative part of the Hamiltonian, we have to 
first rotate % nsi matrix from its flavor basis to the tilde 
basis by. 
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UI^HnSiU23- 



(13) 



4 



Thus our total Hamiltonian {'H = Hm + Hnsi), in 
the tilde basis can be written as a linear superposition of 
the zeroth order Haniiltonian(7Jo) with its perturbative 
parts in that same basis. After the inclusion of the NSI 
matrix, which has its effects at the sub leading part, we 
now redefine our Hamiltonian in the perturbative limit as 
{Hi — > Hi +'Hnsi)- Since the upper bounds of these NSI 
parameters are quite highfSi], therefore we will consider 
two different cases, one with e™^ ~ and the other with 
e™^ ^ e. Once we write the Hamiltonian in the tilde basis 
in this form, we would then look to evaluate the S-matrix. 
The S matrix in the tilde basis is related to the S matrix 
in the flavor basis by, S{L) = U23S[L)ul^, where S{L) — 



Texp 



-i /p dx'H{x) 



and L is the distance traversed. To 



evaluate S{L) perturbatively, we choose f2{x) as f2{x) 
e^^°^S{x), where ^2{x) obeys the evolution equation. 



i — f2{x) = Hin{x). 
ax 



and Hi is written in the form, 



Hi = e'"'"=Hie~'"°\ 



(14) 



(15) 



From ((HI), we would like to deduce f2{x) perturbatively. 
So the solution of the evolution equation followed by 
f^{x), can be written in terms of the Hi matrices as, 



px px px' 

f2{x) = l + / dx'Hi{x') + (-i)^ / dx'Hi{x') / dx" Hi{x") 
Jo Jo Jq 

+ (-^)3 r dx'Hiix') r dx"Hi{x") r dx"'Hi{x"') + 0{e^). 



(16) 



From our previous definition of ^?(a;), we can now write 
the S-matrix as, 

S{x) = e-'""''n{x). (17) 

The S-matrix in the flavor basis is obtained by rotating 
S in the (2-3) space a.s S = U23SUl^. 

Since the S-matrix changes the flavor of a neutrino 
state after traversing a length L, which is given by the 
expression, 

i^ciL) - Sa^^MO), (18) 

the oscillation probability of the neutrino, changing the 
flavor from a ^ P is given as. 



(19) 



treated as a superposition [il, [sll of pure orthonormal 
flavor states. 



(20) 
(21) 




where, e^^ and e'^p are NSI at source and detector, re- 
spectively and the normalization factors are given by 



K = V'[(l+6-^)(l+6^t)]^^, 



0P 



(22) 



This expression of the oscillation probability takes into 
consideration the standard matter interaction and the 
NSI during propagation only. In section HI we will in- 
troduce the idea of NSI at the source and detector. It 
should be noted that since 9i2 and ^23 are quite large, 
compared to ^13, they are considered to be in the zeroth 
order. 



III. NSI AT SOURCE, DETECTOR AND 
DURING PROPAGATION IN v^, 
OSCILLATION PROBABILITY 

In the presence of NSI at the source and at the detec- 
tor, the neutrino states produced at the detector can be 



For example, e^^ describes a non-standard admixture 
of flavor /3 to the neutrino state which is produced 
in association with a charged lepton of flavor a. This 
means the neutrino source does not produce a pure flavor 
neutrino eigenstate Ii'q), but rather a superposition of 
pure orthonormal flavor states To be consistent 

with the literature the convention that we have chosen 
is in e*^, the flrst index corresponds to the flavor of the 
charged lepton, and the second to that of the neutrino, 
while in e^^, the order is reversed. In general, as we can 
clearly see from the definitions above that the matrices 
(1 -h e") and (1 -I- e'^) are non unitary, i.e. the source and 
the detection states do not require to form a complete 
orthonormal sets of basis vectors in the Hilbert space. 
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Since the coefScients e^^ and e^g both originate from 
the {V — A){V ± A) coupling [sTI to up and down quarks, 
so we have the constraint 



IV. PERTURBATION THEORY BY 
CONSIDERING LARGE NSI PARAMETERS 
DURING PROPAGATION 



^ea — ^ae- (23) 

Thus this condition reduces the number of independent 
parameters and makes the model more predictive. But 
in our paper we have presented the most general case. 

Considering the NSI effects at the source, detector as 
well as during the propagation of neutrinos through mat- 
ter the amplitude of the oscillation becomes 

= \{4\s{L)w:,)\' 

= ^Tjhjdli^ + ''f SiL) (1 + en^]p^\%2i) 

where the S{L) is defined earlier. Considering NSI at 
the source and at the detector of the order of e[5j, we 
can now write the probability expression as a sum of 
the probabilities of different order y/e terms. The total 
oscillation probability would look like 
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In our later analysis we would incorporate the results 
from the reactor neutrino experiments along with the 
long baseline superbeam experiment such as T2K. It 
is notable that the non-standard interaction parameters 
during propagation do not play any substantial role in 
case of the reactor neutrino experiments, due to its very 
short baseline. Furthermore, the NSI parameters present 
both at the source and at the detector of these two dif- 
ferent kinds of neutrino experiments i.e. e^^ and e'^^ are 
considered to be the same. The oscillation probability 
P{i'a ~^ vp) is for a neutrino, rather than an antineutrino. 
However, one can relate the oscillation probabilities for 
antineutrinos to those for neutrinos by 



Pap = Pap{5cp -5cp,A^ -A). 



(26) 



In addition, we also have to replace e*, e'', e™^, with 
their complex conjugates, in order to deduce the oscil- 
lation probability for the antineutrino, if one considers 
non-standard interaction during propagation and at the 
source and detector of the experiment. 

It should be noted that the expression ([^^ is also valid 
in the Minimal Unitarity Violation (MUV) model and 
is very instructive for analyzing the CP violating effects 
in the MUV model in future long baseline experiments 

EllilleSil. 



In our next two sections, we will consider two different 
cases of these NSI parameters and present oscillation 
probability for the channel ^ Ve, as was observed by 
the T2K experiment. In this section, we will consider 
/e ^ 0.18, and in the next section, we will put 



^ap 
^ap 



0.03. 



Since A is of the order of e, thus Unsi is in the pertur- 
bative range . Furthermore, wc also have to transform 
Hnsi from its flavor basis to the tilde basis, e.g. 



UnSI — U23'HnSiU23- 



Our total Hamiltonian now looks like 



"H — -ff + Hi + Hnsi 



(27) 



(28) 



We include this Hnsi in the perturbative part of the 
Hamiltonian and follow the same calculations described 
previously. The order e'^/^ component of the Hamiltonian 
now looks like 
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By computing the S-matrix, which comprises of the stan- 
dard matter interaction, and NSI during propagation of 
the neutrino, we then include the NSI parameters at the 
source and detector, which are of the order of e, as per 
the present bounds on these parameters suggest. Finally 
we write down the oscillation probability for the muon 
neutrino going to electron neutrino up-to second order 
in e. It is noteworthy that due to the non unitarity of 
the non-standard interaction matrices at the source and 
at the detector, the Probability of neutrino oscillation is 
not normalized to unity. So one has to include necessary 
normalization factors as we have done in (1221) . In the 
context of T2K, where we are observing muon neutrino 
oscillation to electron neutrino, these normalization fac- 
tors do not play a very significant role. To be precise 
the effects of these normalization terms are greater than 
O(e^), which we are neglecting. Finally by considering 
all these effects, the oscillation probability in the — >■ z/g 
channel is, 
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where S2xij = siii20ij and C2xji = cos20ij and 
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There are a few salient features of this expression of the 
probability. These are as fohows 

• Considering the baseline length to be zero, we are 
left with the term 
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This term is the manifestation of the non-unitarity 
of the source and detector matrices, more com- 
monly know as the zero distance effect. 

Assuming the standard matter interaction and NSI 
during the propagation, as well as at the source 
and at the detector to be absent, we can obtain the 
expression of the probability, representing the vac- 
uum oscillation probability for a three flavor neu- 
trino scenario correct upto 0{a'^) . Particularly the 
following leading term of vacuum oscillation which 

is lei 
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in our expression for the probability. 
For the ~^ Vf. channel, only e™, and el; 



—r t^e "-iiciiiiit-i, "jiiij ciiiii appears 
NSI parameters during the propagation of the neu- 
trino. The contribution from all the other NSI 
parameters during propagation are very much sup- 
pressed. 



• As mentioned, our expression of the probability is 
of the order of e^, by considering lar ge a ngle 6*13. 
Similar expressions are to be found in 5l| , the au- 
thors in there work considered small ^13. But the 
recent reactor based experiments [2, 3] compelled us 
to consider the regime of large sin ^13 ~ y^. 

V. PERTURBATION THEORY BY 
CONSIDERING SMALL NSI PARAMETERS 
DURING PROPAGATION 

In this section we will concentrate on the idea of small 
non-standard interaction parameters during propagation. 
The standard matter interaction is again considered to 
be of the order of e, and we now consider NSI during 
propagation e™^ ~ e. Following the same argument, as 
done in the previous section, the zeroth order Hamilto- 
nian in the tilde basis remains the same. However the 
perturbative Hamiltonian gets rearranged. In this case 
'H.N SI would be of the order of e^. 

The part of the perturbed Hamiltonian (ITU)) , can 
now be written as. 



can be obtained from eleventh term of (pO)) after Hi{e ) 
considering A — > 0. 

• Since muon is produced at the source therefore we 
observe only e*^, e^^, e*,. NSI parameters, and elec- 
tron is obtained at the detector, thus we observe efg. 
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(34) 



We again follow the same procedure as performed in the 
previous section. We computed the S-matrix, after in- 
cluding the standard matter interaction, and NSI during 
propagation. Then we considered the source and the de- 
tector effect. Thus the probability for muon neutrino 
going to electron neutrino for e™^ ~ e is given by. 
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We would again try to emphasize on some of the interest- 
ing features about this oscihation probabihty expression. 

• As expected, we get back the same expression for 
the near detector effect, which was provided in the 
expression ([5^. 

• Similar to the result of previous section, assuming 
the standard matter interaction and NSI during the 
propagation, as well as at the source and at the de- 
tector to be absent, we can obtain the the expres- 
sion of the probability, representing the vacuum os- 
cillation probability for a three flavor neutrino sce- 
nario correct upto 0(0^) . Particularly from the 
tenth term in (1551) considering A — one can get 
the leading vacuum oscillation probability in ([33]). 



• It is very much interesting to note, that due to the 
choice of the NSI parameters during propagation 
proportional to e, the probability expression up to 
second order in e for this particular channel is de- 
void of any terms containing this kind of NSI. It 
shows that it is very difficult to constrain such small 
NSIs in relatively short baseline neutrino oscillation 
experiments like T2K. 



• It is conspicuous, that the NSI parameters at the 
source and at the detector carry the same flavor 
indices, as in Eg. ([3D)) . 




FIG. 1: Plot of eee — S for real NSI in matter. Excluded 
region (white at 90% and Grey+white region at 66% 
confidence level). Upper(lower) panel corresponds to 
normal(inverted) hierarchy. 



paper giving constraints on NSIs. In our numerical anal- 
ysis we have considered even much higher values of NSI 
which has not been considered in our perturbative ap- 
proach but which are allowed after considering model in- 
dependent constraints [5]. Some of the NSIs like eee, £fj,fn 
err (for propagation) etc. do not appear in our expres- 
sion of P^^^y^ in section IV and V as those have been 
assumed to be very small. In our numerical analysis how- 
ever, we still have obtained some constraints on those 
NSIs because of their presently allowed higher model- 
independent values as discussed later. 
For ultra-relativistic neutrinos, we have 

2 4 

Ek^E+^^, pc^E, ct^x. (36) 

where k = 1,2,3 correspond to mass eigenstates and Ek 
are the eigenenergies. The neutrino energy E is the av- 
erage energy after assuming the three momentum of dif- 
ferent components (1,2 and 3) to be equal. 

Therefore, for the numerical analysis the relevant tran- 
sition evolution equation for the flavor transition is 

ihc-^Si3a{x) = '^'HprjSrja, (37) 



with initial condition S'/3q,(0) = Spa and H is the total 
Hamiltonian comprising of standard matter interaction 
and NSI during propagation. However, here apart from 
NSI in propagation we want to include the source and 
the detector NSI interaction. So, to implement that we 
apply the source and detector NSI matrices to the SajS 
which is already mentioned in (j24p as 

The probability expression for the transition Va — >■ is 

P,j^,. = \Apa.\^. (38) 

Although NSI considered at source, that at the detec- 
tor and that during propagation are in general different, 
however, following 43} we have considered e^^ = e^^. 

The recent T2K result [l| has obtained the constraint 
on (5 — sin^ 29i3 plane at 90% confidence level based on 
the events in — > i^g transition in the baseline of 295 
Km. Furthermore, Daya Bay reactor neutrino experi- 
ment has recently measured the neutrino mixing angle 
9i3 with 5.2 a confidence level for which sin^ 29i3 = 
0.092 ± 0.016 ±0.005(syst). Here, we analyze both these 
constraints considering real NSIs (one at a time) in prop- 
agation. Somewhat conservative bound on all NSIs at 
source and detector have been considered and taken to 
be of about 10"'^. To tune with the experimental result of 
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FIG. 3: Plot of Eer - S for real NSI in niatter. Excluded 
region (white at 90% and Grey+white region at 66% 
confidence level). Upper(lower) panel corresponds to 
normal(inverted) hierarchy. 



T2K and Daya Bay we shall use the same range of prob- 
ability of oscillation as one obtains using the constraints 
on 5 — sin^ 20i3 given by T2K (without considering NSI) 
at different confidence level for normal and inverted hier- 
archies (allowing the variation of 6*13 as in T2K paper). 
After that we shall fix O13 at Daya Bay value and find 
out the allowed ranges in the parameters S and different 
NSIs in matter (one at a time) subject to this constraint 
on the probability of oscillation. For numerical analy- 
sis we use the following values as considered by T2K 
[H : Aml2 = 7.6 x IQ-^eY^, Ara\^ = 2.4 x 10~3eV^ 
sin^ 26*12 = 0.8704, sin^ 26*23 — 1-0, an average earth den- 
sity p = Zflgjcw? and central value of sin^ 26*13 = 0.092 
as obtained from Daya Bay. One may note here that 
in presence of NSIs paticularly in matter (which are ex- 
pected to be much larger than those NSIs at source and 
detector) the neutrino mixing parameters considered by 
T2K could change [HjIg^I. However, it is found that only 
when several NSIs are considered simultaneously small 
changes occur in the best-fit values of these parameters. 
As for example, considering solar and KamLAND data 
the change in best-fit value of 6*12 and Am^^ can be seen 
in figure 2 in 66j and considering atmospheric and K2K 
data the change in best-fit value of 6*23 and Am23 '^^^ be 
seen in figure 5 in [63| after considering NSIs like Cee, Cer 
and Ct-t simultaneously. However, in our analysis, we 
have obtained constraint on S— NSI plane by consider- 
ing one of the NSIs at a time for which the changes in 
these mixing parameters are expected to be small and in 




FIG. 4: Plot of e^^ - 6 for real NSI in matter. Excluded 
region (white at 90% and Grey-|-white region at 66% 
confidence level). Upper(lower) panel corresponds to 
normal(inverted) hierarchy. 



this simple analysis we use the values of mixing parame- 
ters as considered by T2K. For a very rigorous analysis 
in obtaining these mixing parameters in presence of NSIs, 
one is required to fit the solar, KamLAND, atmospheric 
and K2K data simultaneously by considering NSIs one at 
a time or all together in the general three fiavor neutrino 
mixing scenario which has not been done so far to the 
best of our knowledge. 

In all the plots of NSI versus S the Eq/j corresponds to 
e™^ which are NSIs in matter during propagation. Dark 
shaded regions correspond to allowed region. T2K con- 
straint on (5— sin^ 26*13 plane at 66% confidence level (C.L) 
together with Daya Bay result on ^13 correspond to ex- 
cluded region (white-|- grey) and only white excluded re- 
gion correspond to the same T2K constraint at 90 % C.L. 
We have done the analysis on NSI constraints keeping in 
view the allowed range of NSIs for earth like matters as 
mentioned in j5|] and have considered those for real values. 
Upper (lower) panel in each plot correspond to normal 
(inverted) hierarchies. Out of various NSIs the significant 
constraints are obtained particularly for e™ and for 
both the hierarchies of neutrino masses. Particularly for 
e^^ and e^,- no constraint can be obtained for normal 
hierarchy. 

In Fig. [T]it is seen that in the upper panel for normal 
hierarchy the constraints on eee can be obtained for neg- 
ative values only for certain values of S corresponding to 
T2K's 66% confidence level result whereas for inverted 
hierarchy the constraint is mainly on positive eee- How- 
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FIG. 5: Plot of e^T- — 6 for real NSI in niatter. Excluded 
region (white at 90% and Grey+white region at 66% 
confidence level). Upper(lower) panel corresponds to 
normal(inverted) hierarchy. 



ever, for inverted hierarchy there is also excluded white 
region corresponding to T2K's 90% confidence level re- 
sult. In Fig. [2] the constraint on is found mainly for 
positive values for both the hierarchies corresponding to 
T2K's 66% confidence level result only. In Fig. [3]eer is 
significantly constrained. For normal hierarchy the neg- 
ative (positive) value could be constrained up to about 
-0.2(0.6) for certain values of 5 corresponding to T2K's 
66% confidence level result. For inverted hierarchy such 
constraints are even more stringent and for certain values 
of 5 all negative values could be excluded. In Fig. |4] no 
constraint is obtained for e^^ for normal hierachy. How- 
ever, for inverted hierarchy around 5 = 7r/2 all values are 
excluded corresponding to T2K's 66% confidence level 
result. In Fig. [5] for normal hierarchy no constraint is ob- 
tained on e^T-. For inverted hierarchy all negative values 
are excluded. In Fig. [B] there are stringent constraints on 
Stt particularly for positive values for normal hierarchy 
and negative values for inverted hierarchy corresponding 
to T2K's result at both 66% and 90% confidence level. 
Significant part of negative (positive) values of err are 
also excluded for normal(inverted) hierarchy correspond- 
ing to 66% confidence level. Interestingly corresponding 
to T2K's 66% confidence level result one may obtain some 
constraints on Et-t independent of d from Fig. [51 




FIG. 6: Plot of — 5 for real NSI in matter. Excluded 
region (white at 90% and Grey-|-white region at 66% 
confidence level). Upper (lower) panel corresponds to 
normal(inverted) hierarchy. 



VII. CONCLUSION 

In this work, using perturbation theory, we obtain the 
probability of oscillation P,y^^^^ (suitable for relatively 
short baseline of T2K and for large 613 as evident from 

Daya Bay experiment) up to order a^{a = ^^i'- ) by 
considering NSIs at the source, detector as well as dur- 
ing propagation of neutrinos through matter. We have 
kept the standard matter interaction part in perturbed 
Hamiltonian which is appropriate for the baseline consid- 
ered by T2K. In addition, we have considered two cases, 
namely, e^^ - e - 0.03 and e^^ ~ ^/e ~ 0.18 -the latter 
corresponding to slightly larger NSI. In the expression of 
oscillation pobability one can see that a flavor transition 
takes place at the source even before the propagation of 
neutrinos due to NSI at source and detector - which is the 
so-called zero distance effect. However, due to stringent 
constraints on these NSI parameters Q we have assumed 
all of them of about 10"'^ in our numerical analysis. Al- 
though there are good model dependent bounds on NSI 
in matter (earth- like) , these are not so strong if one likes 
to constrain them in a model independent way [5|. In 
our numerical analysis, we have obtained constraints on 
various NSIs in matter in a model independent way from 
neutrino oscillation experiments. Nevertheless, one may 
note as mentioned at the end of section V that it is diffi- 
cult to constrain very small NSIs in the relatively short 
baseline oscillation experiment like T2K. Recent Daya 
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Bay result on mixing angle 0i3 has helped us to give 
bound on NSI depending on only one so far unknown pa- 
rameter S in neutrino mixing matrix. Once one knows 
this phase from some short baseline neutrino oscillation 
experiments one may expect the better understanding 
about the possible strength of NSI. Depending on S value 
significant constraint on e™ and could be possible for 
both normal and inverted neutrino mass hierarchies. One 
may note here that in section IV has not appeared in 
the expression of oscillation probability but still we have 
obtained significant constraint on it because of its very 
high presently allowed model independent values Q . Our 
studies indicate that while finding neutrino oscillation pa- 
rameters it might be important to search for any possible 
evidence of NSIs even in the relatively shorter baseline 
neutrino oscillation experiments although longer ones are 
in general preferred. In the coming years the precision 
measurement of neutrino oscillation parameters and the 
NSI parameters in neutrino oscillation experiments could 
be challenging and could even show the evidence of NSIs. 
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